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A new method to construct families 
of complex Hadamard matrices in even dimensions 

D . Goyenechc 

Departamento de Fisica, Universidad de Concepcion, Casilla 160-C, Concepcion, Chile 
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We present a new method for constructing affine families of complex Hadamard 
matrices in every even dimension. This method has an intersection with the Dita 
construction and it generalizes the Szollosi's method. We reproduce well-known 
results, extend some families and present new families of complex Hadamard matrices 
in even dimensions. In particular, found more than 13 millon inequivalent affine 
families of complex Hadamard matrices in dimension 32. We also find analytical 
restrictions for any set of four mutually unbiased bases existing in dimension six. 

Keywords: Complex Hadamard matrices, Affine families. Mutually unbiased bases. 

I. INTRODUCTION 

In recent years, the knowledge about complex Hadamard matrices has exponentially 
increased and several applications to quantum information theory have appeared. They 
are useful to construct bases of unitary operators, bases of maximally entangled states 



and unitary jiepolarisers Complex Hadamard matrices allows to solve the Mean King 



Problem [2H4| , to construct error correcting codes 



es [5 



to study spectral sets and Fuglede's conjecture 
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to find quantum designs 



6| and also 



10| . Furthermore, they are also use 
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ul in 



mathematics for constructing some *-subalgebras in finite von Neumann algebras 
analyzing bi-unimodular sequences, finding cyclic ra-roots 1^ and equiangular lines 

The existence of complex Hadamard matrices in every dimension is assured by 
the Fourier matrices. However, a complete characterization of inequivalent complex 



* Electronic address: |dgoyeneche@cefop.udec.cl| 



Hadamard matrices is known up to dimension five [l^ and the complexity of tlie problem 
suddenly increases in dimension six. This complexity is not only due to the fact that six 
is not a prime power number, like what happened in the mutually unbiased bases prob- 



lem 18|. The problem of determining the complete set of complex Hadamard matrices 
is still open even in lower prime dimensions: it is open in dimension seven, where a one- 
parametric family and a few number of single complex Hadamard matrices are known. 
Also, we do not know if a continuous family of complex Hadamard matrices exists in 
dimension eleven. A complete understanding of complex Hadamard matrices could help 
us solve the Hadamard conjecture and the mutually unbiased bases problem in non-prime 
power dimensions. 

In this work, we present a new method to construct affine families of complex Hadamard 
matrices. This method allows us to find families stemming from a particular subset 
of complex Hadamard matrices existing in even dimensions. This subset includes the 
Fourier and real Hadamard matrices. This work is organized as follows: In Section II, we 
briefly introduce complex Hadamard matrices. In Section HI, we present our method to 
generate affine families in even dimensions. We prove that inequivalent families of complex 
Hadamard matrices stem from inequivalent real Hadamard matrices when our method is 
used. We also find a restriction for any set of four MU bases existing in dimension six. 
In Section IV, we compare our main results with existing constructions: we prove our 
method generalizes the Szollosi method and intersect Di^a construction. In Section V, 
we exemplify our method by constructing families stemming from the Fourier matrices in 
lower even dimension. In Section VI, we construct families of complex Hadamard matrices 
stemming from real Hadamard matrices. We also extend a known family in dimension 
eight and two families in dimension twelve. Finally, in Section VII we summarize and 
conclude. 



II. COMPLEX HADAMARD MATRICES 



In this section, we briefiy resume basic properties of complex Hadamard matrices. 
More detailed explanations can be found in the book of K. Horadam 19|] or in the self 
contained paper of W. Tadej and K. Zyczkowski 20]. A square matrix H of size d is 
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called a complex Hadamard matrix if its entries are unimodular complex numbers and it 
has orthogonal columns. The Fourier matrix defined by its entries 



27ri 



{FaU = e^^^ (1) 

where . = and A: = 0, . . . , ci - 1, is a privileged example because it represents the 
only construction existing in every dimension d [14]. On the other hand, real Hadamard 
matrices, that is, complex Hadamard matrices having real entries, can only exist in di- 
mension of the form d = 4k, where A; = 1/2 or is an integer positive number. The 
Hadamard conjecture states that real Hadamard matrices exist in all such dimensions and 
it represents one of the most important open problems in Combinatorics. Currently, the 
smaller order where a real Hadamard matrix is still unknown is = 4 x 167 = 668. 

We say that two complex Hadamard matrices are equivalent (ifi ~if2) if two diagonal 
matrices Di, D2 and two permutation matrices Pi, P2 exists, such that 

H2 = D^P^H^P^D^. (2) 

In dimensions two, three and five every complex Hadamard matrix is equivalent to the 
Fourier matrix and in dimension four a uniparametric family stems from the Fourier 
matrix. This is the complete characterization of complex Hadamard matrices in d < 5 
jl^ . In dimensions higher than five a complete classification is still an open problem. In 
the case of dimension six, there are several families of inequivalent coniplex Hadamard 



matrices, a three dimensional family being the biggest explicitly known [21 1. 

A complex Hadamard matrix is dephased if every entry of the first row and every 
entry of the first column are equal to the unity. Given a complex Hadamard matrix Hi 
it is possible to obtain a matrix H2 ~ Hi such that H2 is written in dephased form and, 
conversely, if Hi and H2 have the same dephased form then H2 ^ Hi. However, dephased 
form is not unique and we cannot use it to characterize inequivalent complex Hadamard 
matrices. 

In the cases of = 4 and d > 5 there exist arbitrarily close inequivalent complex 
Hadamard matrices; In such cases, we can define a continuous set of them. A construction 



of this kind of complex Hadamard matrices is called a family. A family is affine 



20| if there 



exists a set H{TZ) stemming from a dephased complex Hadamard matrix H , associated 
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with a subspace TZ of the real space of d x d matrices with zeros in the first row and 
column such that 

H{n) = {Ho exp{iR) -.Ren}. (3) 

Here, R determines an m-dimensional space with basis Ri . . . , Rm- We characterize the 
family with the notation H{^), where ^ is an m-dimensional real vector. That is 

HiO = HiRiO) = H o eMmO), (4) 

where R{C) = Yl^i^i-^i- The symbol o denotes the Hadamard product 

{H^oH,),, = iH,),,{H2).„ (5) 

while Exp denotes the entrywise exponential function 

(Exp(/f)),, =exp(/f,,). (6) 

We say an affine family is maximal if it is not contained in any larger affine family 
H{R') stemming from H, where R C R'. If a family is not affine we say it is non-affine. 
Sometimes, a complex Hadamard matrix does not belong to a family. In this case we say 
it is isolated. For example, the _spectral matrix 5*6 and every Fourier matrix defined in 
prime dimensions are isolated |22|. The transpose of a family is still a family. We say 
that two families Hi{^) and H2{v) are cognate if the families H2{i') and H\{^), where t 
denotes matrix transposition, are equivalent. If a family and its transpose determine the 
same family we say it is self-cognate. For example, the four dimensional Fourier family 
f'^'' is self-cognate. 

The problem of determining the maximal family stemming from a complex Hadamard 
matrix is open in dimension higher than five. Even the problem to find the dimension of 
the maximal family is still open. The best approach in order to find this number is the 



defect of a complex Hadamard matrix 
the linear system 



22l | , that is, the dimension of the solution space of 



I 



i?,„o = 0, {0,...,rf-l}, (7) 

Eti H,^kH*{R,^k - Rj,k) = 0, < ^ < J < d - 1, 
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where R G M.'^'^ is a variable matrix. The defect d{H) is an upper bound of the dimension 
of the maximal family stemming from H. For example, d{H) = implies that H is an 
isolated matrix but the reciprocal implication is not valid. The defect has been analytically 
obtained for the Fourier matrices in every dimension 2^ and it has been proven that 
d{Fd) = when d is prime. Therefore, the Fourier matrices in prime dimensions are 
isolated. 



III. CONSTRUCTION OF AFFINE FAMILIES 

A family of d x d complex Hadamard matrices define a continuous set of orthogonal 
bases in when the parameters of the family are smoothly changed. These bases are 
given by the columns of such matrices and they rotates in a very special way, because 
orthogonality is preserved and also every entry of every column is restricted to be a 
unimodular complex number. In this section, we deal with a very particular way of 
rotating columns. The main idea of this work comes from the following question: 

Can a continuous pair of columns of a family of complex Hadamard matrices 
belongs to a two dimensional subspace of ? 

As we will show, this question has only positive answer in even dimensions. Let {4>k} 
be an orthogonal base such that every vector (p^ defines a column of a d x d complex 
Hadamard matrix H. Our goal consists in finding two continuous vectors (paiO ^^'^ 
such that: 

C.l They are a linear combination of 0o ^^^^^ 0i- 
C.2 They have unimodular complex entries. 
C.3 They are orthogonal. 

C.4 The initial conditions 0a(O) = (j^o o-nd 0f,(O) = (pi hold. 

If these conditions are satisfied for any ^ then the vectors {(t>a{.C)^4'b{.C)^4'2^ ■ ■ ■ I'Pd-i} 
define the columns of a family of inequivalent Hadamard matrices H{C,), for any value of 



6 



^ G [0,27r). Let {(fk} be the standard (canonical) base and let 0o and 0i be the first two 
columns of H. Without loosing of generality we can assume that 



and 



^0 

fc=0 



fc=0 

for a given set of unimodular complex numbers {e*"*} such that Ylk=o ^^'^^ ~ ^- Therefore, 
purposing a linear combination of these vectors (C.l) 

</>a(0 = a;(O0o + y(0<^i, 

d-\ 

= + yiOe'"-')^!., (10) 

fc=0 

and imposing unimodular entries in the last equation (C.2) we obtain 

k(e) + y(0e^"1 = l, (11) 

for every k = 0, . . . ,d — 1. These equations are restricted to the condition 

d-i 

J2e''"'=0. (12) 

fc=0 

The coupled system of equations given by Eqs. fllll) and f[T^ have a solution if and only 
if d is an even number and 0i is a real vector. In fact, expanding Eq. (fTT]) we obtain 

+ \y{0\' + 2Re(a;*(02/(Oe^"'=) = 1, (13) 
for every k = 0, . . . ,d — 1. Summing Eqs. f|T3l) and taking into account Eq. f|T2|) we obtain 

l^(OP + b(OP = l, (14) 

for every G [0, 27r) and 

Re(a;*(O2/(Oe*"^)=0, (15) 
for every A; = 0, . . . , ci — 1 and ^ G [0, 27r). A solution of Eq. ffT^ is 

x{0=cos{0 and y{0 = MOe'^^^\ (16) 
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and imposing Eq.f llSp we obtain that 

ei/3{e)e-fc = (17) 

for every ^ 6 [0,27r). Therefore, the only interesting solution found is 

d-1 
k=0 

and an orthogonal vector to (paiO (C.3) is given by 

Mo = ymo-xmi- (19) 

The entries of (pbiO ^lust be unimodular complex numbers (C.2), then 

^im^-i^k = (-1)^^, \fk = o,...,d-i. (20) 

Multiplying Eq.([I7D times Eq.([20D we find that 

e'^(«) = z and e'"^ = (-l)^ (21) 

and consequently 

d-1 

MO = 5^(-l)'e^(-^)'=V;.- (22) 

fc=0 

Here, we can see the condition (C.4) holds. Note that the only difference between the 
pairs {00) 0i} and {4>a{0 > 4'biO} is the exponential term {e^^~^^''^} appearing in Eqs. (fTH]l 
and (122|) . Before formalizing this result in a proposition let us define a useful concept. 

DEFINITION III.l Let Ca and Cb be two columns of a complex Hadamard matrix. 
We say they are an equivalent to real (ER) pair if 

iClUCB)j = ±1, (23) 

for every j = 0, . . . ,d — 1. Here, {CA)j and {Cb)] are the jth entries of Ca and Cb, 
respectively, and the asterisk denotes complex conjugation. 

For example, the Fourier matrices defined in even dimensions d have d/2 ER pairs of 
columns. Indeed, the kth column is given by 

d-1 

(F,), = 5^a;'Vz, (24) 

1=0 
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and the d/2 ER pairs are 

' {Fa)u.{Fa),^a\. (25) 



where fc = 0,...,2 — 1- Given that is symmetric it also has d/2 ER pairs of rows. The 
maximal number of ER pairs of columns and rows of a complex Hadamard matrix H are 
denoted by rj^ and rj^., respectively. These numbers coincide for the Fourier matrices but 



they are different in general. For example, the spectral matrix 5*8 [23[ has r^c = 4 and 
rjr = 0. We formalize the results found in Eqs. (|T8|) and (122|) in the following proposition, 
which represents the main result of this paper. 

PROPOSITION III.l Let H be a complex Hadamard matrix defined in an even di- 
mension d > 2. If H has < m < d/2 ER pairs then H belongs to a m- dimensional 
family of complex Hadamard matrices. 

li m = d/2 a parameter can be missed after dephasing and we do not have a. d/2 dimen- 
sional family. The following corollary clearly emerges from this proposition 

COROLLARY III.l A complex Hadamard matrix having an ER pair of columns or 
rows is not isolated. 

The proof of this corollary is trivial, but our intention here is to remark that, sometimes, 
two columns of H contain enough information to affirm that H is not isolated. Let us 
define a particularly interesting case of ER pairs. 

DEFINITION III.2 Let {C^Ci} and {03,0^} be two ER pairs of columns. We say 
they are aligned if 

iCl)kiC2)k = {C;)k{C,)k, (26) 

for every k = 0, . . . ,d — 1. The maximal number of aligned pairs of columns and rows is 
called rjc and rjr, respectively. 

Let us show that the existence of ER pairs is invariant under equivalence. 

PROPOSITION III. 2 Let H and H be two equivalent dephased complex Hadamard 
matrices and let {Ci,C2} be an ER pair of columns of H. Then, the corresponding pair 
{Ci, C2} is an ER pair of columns of H. 



9 



Proof: Let {Ci,C2} an ER pair of columns of H. That is, 



d-1 



Ci = J2(^M, (27) 



fe=0 



and 



Therefore, we have 



d-1 



C2 = J2{-l)W\^k). {2t 



k=0 



{Ci)kiC2)k = (-1)'. (29) 

Let if be a complex Hadamard matrix equivalent to H. Therefore, there exists unimodular 
complex numbers c and bk, and an injective function f : ^ I'd such that 

d-1 

Ci = ^cbkaf^k)\'fk), (30) 

fc=0 

and 

d-1 

C2 = Y, Sign[(-1)^(°)] c6fe i-iy^''>af^k) \^k). (31) 

fc=0 

The function f{k) and the numbers are related to a permutation operator P and a 
diagonal unitary operator D applied to H, respectively. In order to dephase Ci and C2 
we consider 

c=(6oa/(o))*. (32) 

From Eqs. (l30i) and (l3Ti) we obtain 

(^?i)UC^2)fc = Sign[(-l)^(°)](-l)^«. (33) 

Therefore, Ci and C2 are an ER pair of columns. □ 
Let us present the main consequence of this proposition. 

COROLLARY III. 2 Let H be a dephased complex Hadamard matrix. Then, the num- 
bers ric,rir,rjc and rjf are invariant under equivalence. 

Proof: The numbers t]c and rjr are trivially invariant from the proposition, whereas rjc 
is invariant due to the function / appearing in Eq. fl33|) is the same for every ER pair of 
columns; analogously for rif. □ 
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In the case of the Fourier matrices we have 

if d is even, 

(34) 

if d is odd, 

and the same values for rjr, rjc and rjf. However, for real Hadamard matrices we know that 
Vc = Vr = d/2, but it is not possible to guess the functions rjc and rjf a priori. Let us 
present some properties of ER pairs in dimension six. 

PROPOSITION III. 3 A 6 X 6 complex Hadamard matrix H belongs to the Fourier 
family if and only if r]c ^ or rjj. ^ 0. 

Proof: Suppose that H is written in dephased form and t]c 7^ 0. Therefore, there exists 
H H such that one of its columns isC = (1,-1,1,— 1,1,— 1). Taking into account that 
every entry of H is unimodular and every pair of columns is orthogonal we only obtain 
two free parameters, which generate the Fourier family . Analogously for rj^ ^ 0. □ 
Another way to interpret this proposition is by saying that in dimension six we have 

r]c ^ ^ r]c = r]f = 3, (35) 

and 

r]r ^ ^ r]c = T]f = 3. (36) 

(2) 

It has been proven that from pairs of the form {I, Fg (0^6)} it is not possible to find 



more than triplets of MU bases for any parameters {a,b) 2J]. Therefore, an interesting 
consequence emerges 

COROLLARY III. 3 Suppose that four MU bases can be constructed in dimension six, 
namely {I, Hi, H2, H^} ■ Then, rj^ = rjr = for Hi, H2 and if 3. 

This corollary means that if such matrices Hi, H2, H^ exist, each of them is inequivalent 
to a complex Hadamard matrix having two real columns. 

Let us consider ER pairs in dimension four, where the most general ER pair of columns 
can be written in the form 

Ci = (l,e^^e*^e^^), (37) 
C2 = (l,-e^^e^^-e^^), (38) 
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for any numbers a,b,c E [0, 27r). The — I's can appear in other entries but the above case 
is general up to equivalence. The inner product between these columns is given by 

(^1,^2) = 1-1 + 1- 1 = 0. (39) 

As we can see, the contribution of every pair of entries to the inner product is ±1; this 
is an exclusive property of ER pairs. From Eq. (!39|) . we can show that our construction 
only works in even dimensions; Otherwise, we do not have orthogonal columns. Using 
Eqs. (ITSil and ([22D we construct the continuous pair of ER vectors 

= (e'«,e*("-«),e'(''+«\e'(^-«)), (40) 
C^iO = (e'«,-e'('^-«),e'('+«),-e'(^+«)). (41) 

Dephasing and considering ^ — ^/2 we find that 

Ci(e) = (l,e'('^+«),e'^e'('=+«)), (42) 
C2(0 = (l,-e'("+«),e'',-e'(^+«)). (43) 

Note that ^ only appears in entries such that {Ci)*{C2)j = — 1. This is general easy rule 
to construct families of complex Hadamard matrices from our method. Let us present a 
proposition regarding families stemming from real Hadamard matrices. 

PROPOSITION III. 4 Let Hi and H2 be two inequivalent real Hadamard matrices. 
Then, the families stemming from them using our method are inequivalent. 

Proof: It is enough to show that H2 does not belong to the family stemming from Hi. 
Given the family Hi {C, ) every entry of ^ should be or vr in order to obtain a real 
Hadamard matrix. In this case, it is easy to show that every ER pair of vectors does 
not essentially change. That is, at most we obtain a permutation of rows or columns. 
Therefore, H2 does not belong to ifi(^). □ 
Despite of Proposition IIII.4t Hi ~ H2 does not imply that the families stemming from 



these matrices are equivalent, as we can see in the Dita family D[^j. [25] . 



COROLLARY III. 4 In dimension 32 we can generate 13,710,027 inequivalent affine 
families of complex Hadamard matrices. 
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This corollary is consequence of the recent classification of real Hadamard matrices in 



d = 32, where there exists exactly 13,710,027 inequivalent real Hadamard matrices [26 1. 
In the same way, in dimensions 16, 20, 24 and 28 we can construct 5, 3, 60 and 487 
inequivalent affine families, respectively. Let us present some general cases where our 
construction cannot be applied. 

PROPOSITION III. 5 Let H be a dephased d x d complex Hadamard matrix such that 
one of the following conditions hold 

1. The entries of H are real in the main diagonal and non-real in other cases. 

2. The entries of H are roots of the unity with k = 0, . . . , d/2 — 1, where cu = . 

3. The entries of H are powers of roots of the unity u = , where N is an odd 
number. 

Then, our method defined in Proposition MIL 1\ cannot be applied. 

The proof is trivial because in these cases we do not have ER pairs. □ 
For example, we cannot construct a family stemming from the matrices Cg, Dg, -Dio 
and Du using our method (see the BTZ catalog [27^ for explicit expressions of these 
matrices). Also, this proposition tell us that we cannot construct a family stemming from 
the isolated matrix 5*6 and from every isolated matrix recently found by McNulty and 



Weigert 



28|. 



IV. SZOLLOSI AND DlfA METHODS 

In dimension four or higher than five several families of complex Hadamard matrices 
have been found. Some of them are specific constructions and they cannot be extended 
to other dimensions. A method found by Szollosi allows us to find a family of complex 
Hadamard matrices stemming from real Hadamard matrices. This method is stated in 



Lemma 3.4 



29|: 



Lemma 3.4 (Szollosi) Let H be an arbitrary dephased complex Hadamard matrix of 
order d > A. Suppose that H has a pair of columns, say u and v, with the following 
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property: Ui = Vi or Ui + Vi = holds for every i = 0, . . . ,d — l. Then, H admits an affine 
orbit. 

This lemma is complemented by Theorem 3.5, which represents the main result of the 
Szollosi's paper: 

Theorem 3.5 (Szollosi) Let H he a real Hadamard matrix of order N > 12. Then, 
H admits an (y + 1) -parameter affine orbit. 

We have noted that Lemma 3.4 coincides with our Corollary llll.li Note that the main 
result found by Szollosi is Theorem 3.5, which consider real Hadamard matrices. This 
means that Lemma 3.4 has been defined for complex Hadamard matrices but thinking 
about real Hadamard matrices. Coincidentally, our method is a natural generalization 
of the Szollosi's idea to complex Hadamard matrices. The advantage of our method can 
be appreciated even for real Hadamard matrices, because we can consider ER pairs of 
rows and columns simultaneously. Consequently, we can construct families having more 
than d/2 + 1 independent parameters. In the cases of d = A and d = 8 our constructions 
stemming from real Hadamard matrices agree with Szollosi's results, but in dimension 

(7) " 

twelve our method extends the family found using Szollosi's method, as we will show 
in Section IVll 

In order to construct a family of complex Hadamard matrices we can consider the 



30|. This method can be applied to Di^d type complex Hadamard 



Diva's construction 
matrices: 

Di^a type A complex Hadamard matrix H of order d = did2 is called Di0 type if 
there exists complex Hadamard matrices M of order di and Ni, . . . , N^-^ of order d2 such 
that H can be cast in the form 



( miiNi . . . mid^Nd^ ^ 



(44) 



H = 

y ITLd^lNi . . . ITLd^d^Nd^ J 

The matrices A^i, . . . , A''^^ are not necessarily different. 

Di^a construction Let H be a Di^d type complex Hadamard matrix. Then, the 
following affine family stems from H 
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miiiVi mi2-D2A^2 



mid^Dd^Nd- 



\ 



(45) 



X^md^iNi md^2D2N2 ... md^d^Dd^Nd^ j 



where D2, . . . , Dd^ are diagonal unitary matrices, each of them containing d2 — 1 free 
parameters. The total numbers of free parameters is {di — l){d2 — 1) + m + ni + ■ ■ ■ + nd^ , 
where m and rii, . . . , nd^ denote the number of free parameters of M and A^^i, . . . , A^^^, 
respectively. In the case of Ni = N2 = ■ ■ ■ = Nd-^^ = N we have H = M ® N and Dita 
type coincides with Sylvester type matrices [31.] . Let us present a proposition that shows 
the intersection of our method with the Diva's method. 

PROPOSITION IV. 1 Let H he a dephased complex Hadamard matrix such that rj^ = 
rjf = d/2. Then, H is Sylvester type. 

Proof: Suppose that there exists a complex Hadamard matrix such that rj^ = rjf = d/2. 
Therefore, permuting rows and columns of H we obtain the following equivalent matrix 



where A and B are d/2 x d/2 complex Hadamard matrices. In this case, our method is 
reduced to a particular case of the Di^a construction 



In the case of rj^ 7^ d/2 ot rjr 7^ d/2 our method is not a particular case of the 
Dita construction. This is due to every Di^a-type matrix is embedded within a family 
depending on at least {di — l){d2 — 1) free parameters, where d = did2. On the other 
hand, our method do not have this restriction but we only require the existence of an ER 
pair in any even dimension. We remark it is not easy to identify Di^a type matrices, but 
it is straightforward to identify ER pairs even in higher dimensions. 




(46) 




(47) 



□ 
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V. FOURIER FAMILIES IN EVEN DIMENSIONS 

In this section, we construct families of complex Hadamard matrices stemming from 
the Fourier matrices using our method. We analyze even and doubly even dimensions 
separately because different properties appear. We do not present new results in this 
section because the Fourier matrices have rjc = r]f = d/2 and then our method coincides 
with the Di^a construction, as we have noted in proposition II V. 1 [ The first case we analyze 
corresponds to even dimensions. 



A. Fourier families in even dimensions 



(4J 



In dimension four, the Fourier matrix is given by 

/l 1 1 1 \ 

1 i -1 -i 
1-11-1 
\ 1 -z -1 i J 

Here we have rjc = rj^ = 2 but only a one parameter family can be found after dephas- 
ing. Let (6*2)^ and (6*4)^ be the jth entries of the second and the fourth column of F4. 
According to the practical rule defined after Eq. fH2]) . we should add a phase e*^ in every 
entry of Ci and C2 if (C2)j*(C'4)j = —1- Therefore, we obtain the following family 

F«(e) = F4oexp(zi?^a)(0), 



(49) 



where, 













• 










^ 4 


• 


• • • 












and the symbol • means zero and ^ G [0, 27r) 


. This family 



(50) 



n 

with the only one parametric family existing in dimension four [27]. In the same way, we 
construct a 2-parametric family stemming from the Fourier matrix in dimension six by 
considering the ER pairs of columns {02,0^} and {Cs,Cq}. That is. 



F^^\a,b) = F,oexp{tR';>{a,b)), 



(51) 
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and 



where 



Ff (a,6))* = oexp (^z (i?g(a,6))*^ , 



(52) 



R (2){a,b) 




(53) 



y« a b • a b J 

This is the only maximal affine family stemming from the Fourier matrix Fq. It is easy 
to show that the family Fg^-* (a, b) cannot be extended by considering ER pairs of rows. 
This is a general property valid for d = 4 and even dimension of the form d ^ 4k + 2. 

As we have shown, the Fourier matrices contain d/2 ER pairs of columns in every even 
dimension. Therefore, we can construct the following affine families 



,(d/2-l) 



(0 = i^doexp(zi?W2-i)(0), 



(54) 



and 



where 



Fr~^^(0)* = F.oexp (.(i?(^/-^)(0)*) 



(55) 



• when z = or even 
)) . =<; • when i is odd and j = mod {d/2) (56) 
otherwise. 

The above families are not self-cognate. We have noted from the BTZ catalog {27] that in 
the case of ci = 2, 4, 6, 10, 14 our construction agrees with the maximal affine Hadamard 
family stemming from Fd- This motivated us to establish the following conjecture 

CONJECTURE V.l The maximal affine family of complex Hadamard matrices stem- 
ming from the Fourier matrix Fd in dimensions d = 2p (p prime ) is given by Eq. ( [5^ to 
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B. Fourier families in double even dimensions 



In the case of dimensions of the form d = 4k, A; > 1 we can simultaneously apply 
Proposition IIII.ll to rows and columns of the Fourier matrices, increasing the number of 
parameters of the family beyond d/2 — 1. For example, we have obtained the following 
5-parametric self-cognate family stemming from Fg 



F^^\a, b, c, d, e) = Fgo exp{iRpg{a, b, c, d, e)). 



(57) 



where 



Rpsia, b, c, d, e) 




{51 



• a + d e a • a + d e a 

• b • b • b • b 

\» c + d e c • c + d e c / 

The ER pairs here considered are {(^2, Ce}, {C3, Cy}, {C4, Cg} and 
{R2,Rfi},{R3,R7},{RA,R8} or C = R = {2, 6; 3, 7; 4, 8} for short. The ER pair 
{C4, Cg} produces a linearly dependent parameter and it has not been considered in 
Eq.f l58p . The 5-parametric family given in Eq.f l57p coincides with the maximal affine 
family stemming from Fg In the case oi d = 12, the following 9-parametric family 

can be found 



Fi2\a,b,c,d,e,f,g,h,i) = F^ o exp{iRFr,{a,b,c,d,e, f,g,h,i)), 



(59) 
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where 



RF^2{0',b,c,d,e,f,g,h,i) 



( 

• a+fga+hi a 

• b • b • b 

• c+fgc+hi c 

• d • d • d 

• ^ + f 9 ^ + h i e 

• a+fga+hi a 

• b • b • b 

• c + f g c + h i c 

• d • d • d 

y» e + fge + hi e 



\ 

a + f g a + h i a 

b • b • b 

c + f g c + h i c 

d • d • d 

e + f g e + h i e 

a + f g a + h i a 

b • b • b 

c + f g c + h i c 

d • d • d 

e + f g e + h i e j 



(60) 

Here, we considered the ER pairs C = R = {2,8; 3, 9; 4, 10; 5, 11; 6, 12}, where {Cq, C12} 
produces a hnearly dependent parameter. This result coincides with the affine family 
Fi2A 0i which is self-cognate. A 13-parametric (non-affine) family stemming from F12 
has been recently found by Barros and Bengtsson 32[. This family contains the set of 
affine families found by Tadej and Zyczkowski [20] and, consequently, also contain our 
family F^f^ . Finally, we obtained a 13-parametric family stemming from Fiq. That is. 



^16^^ (a, 6, c, d, e, /, g, h, i, j, k, /, m) = Fig o expiiRp,,, (a, 6, c, d, e, /, g, h, i, j, k, l,m)), (61) 
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where Rp^^ (a, b, c, d, e, f, g, h, i, j, k, I, m) is given by 



v 



a + h i a + j k a + I m a 

b • b • b • b 

c + h i c + j k c + 1 m c 

d • d • d • d 

e + h i e + j k e + I m e 

f • f • f • f 

g + h i g + j k g + I m g 



a + h i a + j k a + I m a 

b • b • b • b 

c + h i c + j k c + l m c 

d • d • d • d 

e + h i e + j k e + I m e 

f • f • f • f 

g + h i g + j k g + I m g 



\ 



a + h i a + j k a + I m a 

b • b • b • b 

c + h i c + j k c + 1 m c 

d • d • d • d 

e + h i e + j k e + I m e 

f • f • f • f 

g + h i g + j k g + I m g 



a + h i a + j k a + I m a 

b • b • b • b 

c + h i c + j k c + 1 m c 

d • d • d • d 

e + h i e + j k e + I m e 

f • f • f • f 



(62) 



g + h i g + j k g + I m g J 

This family is self-cognate. The ER pairs considered here are C = R = 
{2, 10; 3, 11; 4, 12; 5, 13; 6, 14; 7, 15} and the ER pair {Cs, Cie} produces a linearly depen- 
dent parameter. We remark that the maximal affine family stemming from the Fourier 



matrix in dimension 16 has 17 parameters 



30|. In general, for d = 4k, k > 1 we 



can construct the following families of complex Hadamard matrices stemming from the 
Fourier matrices 

Fd{a, P)=FdO exp{t{Ri{a) + R20))), 



and 



where 



Fd{a, /3) = Frf o exp (i (R,{a) + i?2(/3)' ^ 



(63) 



(64) 



• when z = or even 

• when i is odd and j = mod {d/2) 
a[j-i] otherwise, 



(65) 



20 



and 

R20) = Rl0). (66) 

— * 

Note that the parameters a, ^ can be dependent and d — 2 is an upper bound for the 
dimension of the family that we can obtain from Fd. Curiously, we found a d — 3 dimen- 
sional family in the cases of o? = 4, 8, 12, 16 and, maybe, this is also be valid for every 
double even dimension. 

Using our construction, we can define a lower bound for the maximal affine family 
stemming from the complex Hadamard matrices in even dimensions. 

PROPOSITION V.l Let H he a complex Hadamard matrix and dmax{H) the dimen- 
sion of the maximal affine family stemming from H. Then, the following lower bound can 
he established 

Vmax{H) < dmax{H). (67) 

where rj^ax is the maximal number of linearly independent parameters of a family con- 
structed using our method. Let us analyze this bound in the case of the Fourier matrices 
in even dimensions. 



d 


^max (^of) 


dmax (^d) 


d 




dmax (-^(i) 


2 








10 


4 


4 


4 


1 


1 


12 


9 


9 


6 


2 


2 


14 


6 


6 


8 


5 


5 


16 


13 


17 



As we can see, up to d = 14 every known affine family stemming from the Fourier matrices 
can be constructed using our method. 

VI. REAL HADAMARD MATRICES 

From real Hadamard matrices and the Fourier matrices we obtain d/2 ER pairs of 
columns and rows. However, the ER pairs are not necessarily aligned for real Hadamard 
matrices, and therefore we have many ways to choose them. This is a remarkable difference 
between real and Fourier Hadamard matrices. As we will show next, in the case of d > 8 
the ER pairs should be intelligently chosen in order to maximize the dimension of the 
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families. In li = 2 we have an isolated real Hadamard matrix and for d = 4 every 
real Hadamard matrix belongs to the one parametric family F^^\^) presented in the 
previous section. So, the first interesting case here corresponds to c? = 8, where every real 
Hadamard matrix is equivalent to 



(l 
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1 


1 


1 


1 
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-1 
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-1 


-1 


1 


-1 


1 


1 


1 


-1 


-1 


-1 


-1 


1 


1 


-1 


1 


1 


-1 


1 


-1 


-1 


1 


-1 


1 


-1 


1 


-1 


1 


-1 


1 


-1 


-1 


1 


1 


-1 


-1 


1 


1 


1 


-1 


-1 


1 


1 


-1 


-1 




-1 


-1 


-1 


-1 


1 


1 





(68) 



In order to maximize a family obtained from ifg, we should choose the number of ER 

pairs of rows such that the ER pairs of columns is maximal. Therefore, wc obtain a family 
having 8 free parameters but only 5 of them are linearly independent after dcphasing. 
Applying our method to rows of and considering 77,^ = 4 we obtain, as a first step, the 
4-parametric family 

Hs{a^ b, c, d) — Hso exp(ii?ifg(a, b, c, d)), 



(69) 



where 



RHsia,b, c, d) 





a 


• 


a 


• 


a 


• 




• 


b 


• 


b 


• 


b 


• 


b 


• 


c 


• 


c 


• 


c 


• 


c 


• 


c 


• 


c 


• 


c 


• 


c 


• 


a 


• 


a 


• 


a 


• 


a 


• 


d 


• 


d 


• 


d 


• 


d 


• 


d 


• 


d 


• 


d 


• 


d 




b 


• 


b 


• 


b 


• 





(70) 
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Therefore, we have 9 inequivalent ways to choose the ER pairs of cohimns. That is, 





Ca = 


{1,3;5,7;2,4;6,8}, 




Cb = 


{1,3; 5, 7; 2, 6; 4, 8}, 




Cc = 


{1,3; 5, 7; 2, 8; 4, 6}, 




Cd = 


{1,5;3,7;2,4;6,8}, 




— 


|i, 0, -3, i , Z, D, 4, S|, 




Cf = 


{1,5; 3, 7; 2, 8; 4, 6}, 




Cg = 


{1,7;3,5;2,4;6,8}, 




Ch = 


{1,7; 3, 5; 2, 6; 4, 8}, 




Ci = 


{1,7; 3, 5; 2, 8; 4, 6}. 


Let us exphcitly construct the 


case C 


A- That is. 




c, d, e) 


= Hso exp{iRHg^ia, b, c, d, e)) 


where 








• 


• • • • • 



Rhsa ^) c, d, e 



a+b+d • a+b+d d a+b+d d a+b+d 
a+c+e • a+c+e • a+c • a+c 
a+c+e • a+c+e • a+c • a+c 



V 



(71) 
(72) 
(73) 
(74) 
(75) 
(76) 
(77) 
(78) 
(79) 

(80) 



a+d+e • a+d+e d a+d d a+d 
a+d+e • a+d+e d a+d d a+d 
a+b+d • a+b+d d a+b+d d a+b+dj 

We can considerer 7!! = 7x5x3xl = 105 different choices for the ER pairs of rows 
and the above 9 cases represent only one of these choices. In principle, we do not know 
how many of these families are inequivalent. Using Pita construction, it was found 9 
inequivalent 5-parametric families stemming from Hg 33|, which are a particular subset 
of our solutions. This is very easy to show from Proposition IIV.II 

In Section HVl we mentioned our method generalizes the Szollosi's method even for real 
Hadamard matrices. Let us present an example where this is clearly showed. From the 
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real Hadamard matrix 



H 
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parametric family from considering 


the 


ER 


pairs 



(82) 



{1, 6; 2, 3} and R = {1, 5; 2, 11; 3, 7; 4, 8; 6, 12; 9, 10}. That is, 

H[l\a, b, c, d, e, f, g, h) = H12 o ex.p{iRH^^{a, b, c, d, e, f, g, h)), 

where Rhi2 ('^j d, e, /, g, h) is given by 
/ 



(83) 
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(84) 
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We have checked that all the parameters are linearly independent. Considering different 

combinations of ER pairs of rows and columns we can generate 12!! = 46080 families, but 

(s) 

we do not know how many of them are inequivalent. Our family is a new result that 
extends the family if^^ found using the SzoUosi's method 29|. 

We have noted that our method can be applied to other complex Hadamard matrices 



apart from the real and the Fourier matrices. For example, from HIq (25| we found the 



following 5-parametric family 



i7;[o\a, b, c, d, e) 



h[q^ o exp(zi? {u;)(a, b, c, d,e)), 



^5) 



where 
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(87) 



In this case, we have t]c = 5 and rjr = 5, but they cannot be simultaneously considered 
in order to obtain more free parameters. It was proven that a 7-parametric family stems 
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from h[q^ [25], and Eq.f lHSj) represents a subset of this family. We have also extended 
the 7 parametric family -0^22 |2a| found by using the Dita construction. This extension 



is obtained by considering the ER pair given by the first and the last row in every sub- 
family of -D^2S (straightforwardly obtained from inspection of the family 34 1)- We have 



proven that the new parameter is linearly independent to the rest of the parameters after 

(7) 

dephasing. We have omitted details here to abbreviate, but the twenty subfamilies -Di2A 

(7) 

to -Di2ij generalized in the same straightforward way. As a last result, we showed that 
the single matrix Du presented by Dita |25| 
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belongs to the intersection of our extension of -Di2S5 namely -Di2e- That is. 



R 



Di2 G fl D 



(8) 

i2r- 



59) 



T=A 



Indeed, every subfamily of -Di2S stems from H12 (defined in Eq. fl82|) ) or from a real 
Hadamard matrix equivalent to ifi2, and all these equivalent matrices have the same first 
and last row as H12. In order to obtain D12 we start multiplying from the second to the 
seventh column of H12 by i times and we introduce a parameter ^ by applying our method 
to the ER pair {_Ri, Ri2}- In the case of ^ = — 7r/2 we obtain D12. Analogously for every 
subfamily -D^^r, F = A, . . . , i?. 
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VII. SUMMARY AND CONCLUSION 

We presented a new method to construct families of complex Hadamard matrices in 
every even dimension by introducing the concept of ER pairs. We have reproduced some 
previously known results: 

(i) Maximal families stemming from the Fourier matrix in d — A, 6, 8, 10, 14. 

(ii) Subfamilies stemming from the Fourier matrix in every even dimension. 

We have extended the families that can be constructed from every real Hadamard matrix 

(Hi) Our method generalizes the Szollosi's method, in the sense that the number of free 
independent parameters obtained from the Szollosi's method in d > 12 can be increased 
using our method. 

Our method has an intersection with the Diva's construction: 

(iv) If a d X d complex Hadamard matrix has d/2 aligned ER pairs then it is Sylvester 
type and therefore Diva's construction coincides with our method. 

We have extended two families in dimension twelve: 

(v) We found 46,080 different ways to generalize the family obtained from Szollosi's 
method. However, we do not know how many of our extensions are inequivalent. 

(vi) The family D^^^t, obtained from Diva's construction was extended to D^^y,- Also, we 
have proven that the single matrix D12 presented by Di^a belongs to every subfamily of 

^(8) 

We have established a connection between the mutually unbiased (MU) bases problem in 
dimension six and the ER pairs: 

(vii) Let {I, i^i, i72, -f^s} be a set of four MU bases existing in dimension six. Then, 
Hi, H2 and H^ do not have ER pairs. 

And finally, from our method we generated inequivalent affine families from inequivalent 
real Hadamard matrices. For example: 

(via) In dimensions 16, 20, 24, 28 and 32 we can construct 5, 3, 60, 487 and more than 
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13 millon inequivalent families, respectively. 

We hope to generalize the concept of ER pairs in order to introduce parameters in 
three or more columns of a complex Hadamard matrix. This extension could allow us to 
construct affine families of complex Hadamard matrices in every dimension d > 3. 
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